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ABSTRACT: We consider the behaviour of fermions in the background of 
instanton-antiinstanton type configurations. Several different physics problems, from 
the high energy electroweak interactions to the study of vacuum structure of QCD 
and of large orders of perturbation theory are related to this problem. The spec- 
trum of the Dirac operator in such a background is studied in detail. We present an 
approximation for the fermion correlation function when the instanton-antiinstanton 
separation (i?) is large compared to their sizes (p). The situation when the instanton- 
antiinstanton overlap and melt, is studied through the behaviour of the Chern Simons 
number as a function of R/ p and X4. Applying our results to widely discussed cases 
of fermion-number violation in the electroweak theory, we conclude that there are 
no theoretical basis for expecting anomalous cross sections to become observable at 
energies in 10 TeV region. 
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1. Introduction. 



The purpose of the present paper is to study the behaviour of chiral fermions in 
the background of instanton-anti-instanton type. In particular we wish to understand 
how unitarity works out in the presence of topologically nontrivial effects such as 
instantons, treated within the semiclassical approximation. 

The interest in this problem arises at least from three different sources. 

First, there does not seem to be a universal consensus as yet among the physicists 
whether the fermion-number violation becomes strong in the TeV region scatterings, 
induced by the instanton or sphaleron in the SU(2) electroweak interactions [1-30]. 
Both direct calculations of the cross section by the instanton method as in the works 
of 1^ , and another approach |TB|, which makes use of the optical theorem and 
the so-called valley method, encounter various technical difficulties. The opinion in 
favour of the anomalous processes becoming observable, has been expressed based on 
some toy-model calculations []15| , or on some general considerations |^ . In the valley 
approach the question of unitarity and the behaviour of the fermion propagators in 
the background of instanton-anti-instanton background are quite central. The results 
of our investigation show that there are no theoretical basis in the claims made in 
the literature that the anomalous cross sections become observable in the TeV region 
scattering in the standard electroweak model.Q 

Secondly our study is closely related to that of the vacuum structure in QCD, es- 
pecially in connection with the so-called instanton liquid model although in the 
latter one needs the knowledge of fermion propagators in the multi-instanton back- 
ground. A detailed study of the fermion propagation in the instanton-anti-instanton 
background should help understanding of a more complex situation in the multi- 
instanton media. 

Finally the contribution of instanton-anti-instanton to the correlation functions 
has a deep connection with the problem of large order perturbations both in QCD 
and in the electroweak theory. Many related issues (Borel summability, determination 
of the large order coefficients for R, the role of renormalons , etc.) are receiving a 



^ There are also arguments |l^, |T^, |T^, |3^, which rely on the unitarity constraints in multiple 

gauge-boson productions, that such processes remain necessarily suppressed by a finite fraction of 
the 't Hooft factor, ^ exp(— 47r/a). As far as the authors know these arguments have never been 
rebuked, but at the moment it seems to be difficult to make them more quantitative and rigorous. 
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renewed interest in the literature |3^, |38 



To be definite we consider an SU{2) gauge theory with Np left-handed fermion 
doublets. 

The main problem regarding unitarity and chiral anomaly can be formulated as 
follows. 

The optical theorem states that the cross section, 

1 + 2 — >X, (1.1) 

summed over X, is equal, apart from a kinematical factor, to the imaginary part of 
the forward elastic amplitude, 

1 + 2 — ^1 + 2. (1.2) 

Now consider a particular class of processes ( |1 . 1| ) induced by an instanton ("anoma- 
lous processes"), with the change of the fermion number, 

Af = f, + f2-fx = NF. (1.3) 

Sum over the final states satisfying ( |1.3|) should give a contribution to the imaginary 
part of the elastic amplitude. 

For an (anti-) instanton background, which is relevant for the calculation of the 
production process ( |1 . 1|) , each right (left) handed fermion field has a zero mode. The 
standard functional integration over fermions yields a product of these zero modes; 
by going to momentum space and by applying the LSZ amputation one finds the 
S-matrix elements consistent with ( [1.31 ). 

The corresponding contribution in the elastic amplitude ( |1.2| ) must arise from a 
sort of instanton-antiinstanton {i — a) background, topologically (globally) equivalent 
to the trivial, perturbative vacuum. One expects however no fermion zero modes to 
exist in such a background (see Section 2 for more details). How can one compute 
the "anomalous" part of the elastic amplitude then? 

^The work of Ref. (see also Ref. [^) goes quite some way in proving unitarity in the presence 
of instanton effects, especially as regards the final states (i.e., states summed over in the unitarity 
relations). However, fermions are not considered there: as a result no discussion on subtleties related 
to chiral anomaly is found in [pl]| . In any case, the central issue related with the initial (or external) 
states has not been addressed before. 
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On a general physical ground one expects that at large i — a separation (i?) the 
standard fermion zero modes in the antiinstanton or the instanton background should 
play an important role. For instance the two point function, 



I{x,y) = J Vi/jVip 'ilj{x)'ijj{y) exp — J d xiijjDi/j, (1.4) 

is expected to behave approximately as: 

Iix,y)c^4%)C%r (1.5) 

where r]Q\x) {Co\y)) is the left-handed (right-handed) zero mode in the antiinstanton 
(instanton) background. Such a behaviour was simply assumed in the existing liter- 
ature so far. We wish however to compute I{x,y), prove that Eq.( |1.5|) indeed holds 
approximately at large R, and to calculate the corrections. To do this requires a 
systematic way of calculation. 

To make the problem well-defined, we consider a particular class of i — a type 



configurations - the so called valley, or streamline, trajectory |2^, \L5 



J^ivalley) ^ ^(a) ^ ^(i) ^ ^{int.) ^ ^^^g) 



9 {x - XiYiix - Xi)^ + p^) 

A (int.) _ rr n _ A \ r ~ ~^ y)u _ (x — Xj)^ -, 

' " g^'"" '''''^\x-x, + yY (x-x,)2J' 

where 

y = -R/{z - 1); z = {R^ + 2p' + ^ R^ + Apm^)/2p^; R>' = {x, - xj^ (1.7) 

There are several reasons for such a choice. First, the classical field Eq.( p..6|) is 
known, at least at large R, to have the correct form of interaction between the instan- 
ton and antiinstanton (dipole-dipole). Secondly, Eq.( |1.6|) interpolates two solutions 
of the Euclidean field equations, = A^^'> + A^^'> (at R = \xa — Xi\ = oo) and 
(gauge-equivalent of) = (at i? = ), so that the phenomenon of i — a "melting" 
at -R — > can be studied quantitatively. Thirdly, they are solutions of the valley 
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equation |jT5|. f\ Importance of the non-Gaussian integrations along such an almost 
fiat valley in the field configuration space, was emphasised first in Ref . [p9|, |30| in a 
general context of quantum mechanics and QCD. Furthermore, the valley field of 
Eq.(p..6D satisfies the simple, covariant gauge condition, 

SO that all calculations can be done in a manifestly covariant fashion. 
We shall study the four point function, 

< Ti)i{x)'ilj2{u)ipi{y)'4)2{v) >(^''-"=y) 

P^P^V'i(a:)V'2M^i(y)^2(t;)e-Vz(^=°); (1.9) 

Np 



in the fixed background of Eq.( |1.6| ). Integrations over the collective coordinates such 
as R and p are to be performed afterwards. 

As the functional integral factorises in flavour we must study (suppressing the 
flavour index). 



I{x,y) = J Vijj'Dipip{x)^jj{y) exp— J d xii/jDtjj, (1-10) 

and 

Z = JV'ipV'ipexp-Jd^xiipDip = detD, (1.11) 

where it is assumed that det D is suitably regularised. 

The paper will be organised as follows. First we study in Section 2 (also in 
Appendix A and Appendix C) the spectrum of the Dirac operator in the valley back- 
ground Eq. ( |1.6D in detail. Among others, the absence of the fermion zero modes (and 
■^This is true only in unbroken gauge theories such as QCD. In the presence of Higgs fields 



phenomenological relevance of the valley configuration Eq.(1.6) is not obvious. We concentrate here 



on the problem of unitarity versus anomaly as formulated above, in the (relatively) simple setting 



of Eq.(1.6), which is the starting point also for similar problems in the standard electroweak theory. 

*We believe that, for a well-separated instanton-antiinstanton pair, the valley method does give 
the dominant contribution to the anomalous imaginary part of amplitudes. For < 1, instead, we 
see no reasons to expect such a valley field to be physically distinguishable from generic perturbations 
around = (see Section [4.1]), hence to be of particular importance in the functional integration. 



Nonetheless the use of a concrete and explicit valley field such as Eq.(1.6) is quite adequate for our 
purposes. 



4 



actually of any normalisable modes) is proved. In Section 3 the behaviour of the 
fermion Green functions is studied at large R/ p, leading essentially to the behaviour 
Eg .( 11.51) . In Section 4 we study the situation at small R/ p where the instanton and 
anti-instanton overlap on each other and melt. Although this analysis is somewhat 
indirect, being based on the behaviour of the Chern Simons number for the back- 
ground Eq.( p..6|) as a function of R/ p and of X4, it leads to a semi-quantitative idea 
of when the instanton-anti-instanton pair actually melt. 



2. Spectrum of the Dirac Operator in the Valley Background 

We first want to learn all we can about the Dirac operator V in the valley back- 
ground and in the full euclidean space R^, where: 



V 



I?_ 
V+ 



V+ = iD = io-f^Df,, = iD = ia^D^ 

and = {d — igA)^ is the usual covariant derivative. The discussion will be divided 
in three parts, concerning the essential spectrum 0, the zero modes, and finally the 
positive (normalisable) modes. 

Let us consider the valley field 

Aa = cr^{afj,au - o^ivlf^r — 5- 

g 

= 7 \2 I 2 ~ 7 \^~r^ {v = x-Xi + y) 

[X — Xa) + P [X — XiY + 



which is obtained from Eq. (|1.6|) by the gauge transformation: U = (which 
obviously does not change the spectrum of our operator). The advantage of this 
gauge choice ("clever gauge") is that v4^(x) is now a bounded function Q and that 
A^{x) ~ as X ^ 00, (which implies the uniformity of convergency in the direction 
of X.) We show in Appendix A that these properties allow us to prove the theorem 

o{V) = OessiV) = (-00, +00) : 



^The usual definition of essential spectrum of an operator (see [Q) is aess = '^/o'disc, where 
the discrete spectrum adisc is the set of isolated point of the spectrum that are eigenvalues with 
finite multiplicity. 
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If we define || / ||oo= sup^|/(a::)| (i.e. the minimum upper bound of |/(a;)|), then || A^^ \\oo< 00 
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the essential spectrum is the same as that of the free Dirac operator. 

We now turn our attention to the possible zero modes of the Dirac operator in 
the valley background. Thanks to the property -F^ < oo of our vector potential, 
we can consider (see |Q) the Dirac operator trasported by stereographic mapping in 
the compactified space S"^ (which we call T>), and apply the Atiyah-Singer theorem 



49| to evaluate the index of differential operators in compact manifolds to V. This 



yields the well known result : 

2 

iV+ - iV_ = -— / Tr(FF) 
+ 167r JRi ^ ' 

where N± = dim(kerl?-|-), i.e the number of normalisable (in the sense of measure) 
zero modes of the operator with left (resp. right) helicity. It can be proved that zero 
modes ip of V can be mapped into zero modes ip oiV (and viceversa) and that the 
normalisation condition for in gives the condition 



(i.e. they are in /.^(R'', ^^^^2) )); which does not guarantees the usual normalisation 
condition. Clearly a zero mode of V normalisable in the usual sense can be mapped in 
a normalisable zero mode of T>, so that we have N± < N±, where N± = dim(kerP-i-), 
i.e the number of normalisable zero modes of euclidian Dirac operator with left (resp. 
right) helicity. In the valley background /j^4 Tr{FF) = 0. The question is: is "0=0-0" 
or "0=1-1" or worse ?, i.e, are there zero modes in the valley background ? 

Because the valley is essentially the sum of istanton and anti-istanton background 
(each with one zero mode, but with opposite helicity) one expects that V can have 
at most two zero modes with different helicity. The following assertion will now be 
proved: 

"0 7^ 1 - 1. " □ 

Proof: Let us go back to Lorentz gauge of Eg. (|1.6|): 

= --W^a, - 5^.) F, = fil^F.a'' (2.1) 



F,[x) ^\d, log L[x) L{x) ^ [^''il^^l (x - X, + yr. (2.2 

Z \^X Xi J \ P 



^The claim is valid also for zero modes not normalisable in the usual sense, but that arc in 
L^(R^, -p^ifsy), it suffices to modify the normalisation condition, without other changes. 
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We assume (without less of generality) that Xi = (—-|, 0,0,0) and Xa = (-^,0,0,0), 
so that also the parameter y above has the form y = {y, 0, 0, 0) (see appendix D for 
definition of y) . 

It turns out to be simpler for our purposes to work in a noncovariant formalism: 
^ -dt + ia -V + if ■ F - a At ■ F - Fqg ■ f 
V+^+dt + ia-V -ir-F-aAT-F-Foa-T 

where = {t,r), 

F = rf{r,ty, Fo = Fo(r,t) 

r = V^, (J are spin operators, and r are isospin one. First let us note that V 
commutes with the threedimensional total angular momentum operator J = —if A 
V + c? + r , so that any eventual normalisable eigenmode of V can be chosen to be 
also eigenmode of and J3, with quantum number j, m, as usual. Note that the 
symmetry forces eigemodes to be in 2j + 1-plets, so if we find a zero mode with total 
spin j then there must be 2j more zero eigenmodes. Since we are interested in a 
single zero mode each of and for the reasons mentioned above, we look for a 
solution with j = 0. Let us concentrate on the lefthandcd mode, associated with 
below. 

There are two ways to form a j = state, because the fermion wave function 
is in a doublet representation both of isospin and of spin. Thus, a and r combine 
into a singlet or triplet representation, coupled with a singlet or triplet representation 
respectively of spatial angular momentum, to form a singlet of the total spin. The 
most general function with j = can be written accordingly as: 

Vja = {cr2)jaS{r,t) - i{B(T2)ja ' xT{r,t). 
The action of V_^. on 77 is: 

V+ri ^a2{dtS + 3FoS-x-VT-3T + F-xT) + iaa2 ■ {-VS - 3FS - xdtT + xF^T) . 

Imposing "D+Ty = 0, we get the system of equations: 

dtS + 3FoS-x-WT-3T + F-xT = 0, 

-VS -3FS -xdiT + xF^T = 0. (2.3) 
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Eq.( p.3D can be simplified by defining two new functions T(r,t) and S{r,t): 

T = ^L^f] S = L-IS, 

{L{x) is defined in (|2.2| )), so that Eq. (^l3| ) reads: 

dtS-hdrf =0; 
~drS-hdtf =0; (2.4) 

where 

h = h{r,t) = ^^. 

The consistency condition drdtS = dtdrS following from Eq. (|2.4|) is equivalent to 
V(/iVT) = (here V = {dr, dt)). Furthermore the normalisation condition / |?7p < oo 
imposes separate conditions on 5* and T: 

r+oo r+oo r+oo r+oo T _ 

<fxAT\^ = dt dr 47rr^|Tp = dt dr 4n—T^ < oo; 



OO -/Q J— CO Jo 

//'+00 /'+00 /'+00 r-\-oo _ 

rf1r|5p = / dt dr 47rr2|5p = dt dr Anr^L-^S^ < oo. 
J-oo Jo J-oo Jo 

Note that with our choices of Xi = (— f , 0, 0, 0) and Xa = (f , 0, 0, 0), we have L(r, t) ^ 
almost everywhere on the line r = 0, so that a necessary condition for normalisation 
is T(0,t) = 0. 

We are thus led to the boundary value problem: 

V(WT) = 0, on n; 
f = 0, on 
T G H'^iSl), f G C°(^]), (2.5) 

where f2 = {(r, t) G R^/r > 0} (see Appendix A for definition of H'^). The operator 
+ V/iV is an elliptic differential operator with coefiicents analytic in fi, so that 
we can apply: 

Theorem (Elliptic Regularity) : Let P{x, d) an elliptic differential operator in Q with 
coefficients analytic in Q and let f a function analytic in Q. If u & V{Q) (the space 

®A differential linear operator of order m P{x,d) — X]|a|<m '^"(^)^" called elliptic if 
E|a|=™a«(a;)p" ^Oyx en Vp e R"/{0}; it is called strongly elliptic if E|al=m a"(a;)p" > 
OVxen Vpe R"/{0}. 
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of distributions) is a solution of P{x,d)u = f then also u is analytic in Vl (see 
p.178]). 



T must thus be an analytic function in fi, hence the following theorem (see 
p. 263]) applies: 

Theorem (Strong Maximum Principle) : Let P = aijdidj + bidi be a strongly elliptic 
operator with real continuous coefficients on an open connected region Q C R" then: 
\/u G C^(fi) such that Pu = on i7 and u non constant on Vt, 

inf^M < u{x) < sup^M, Vx G Vl. 

In other words, T cannot have global extrema on f2; using the continuity on H, the 
boundary conditions Eq.( p.5|) and the condition T ^ when r, t ^ oo, it is easy to 
prove that T = 0. 

It follows that S must be a constant, i.e, 

3 

S = const. L(x)~2 . 

The normalisation condition however forces S* = 0. Thus no left handed normalizable 
zero modes exist. Q.E.D. 

Having disposed of the zero modes, one might then ask if the zero mode of P4. in 
the antistanton field, and that of in the istanton background, combine in some way 
in the valley field so as to form a pair of non zero eigenvalue ±A of T) (a sort of quasi 
zero mode, approaching zero as the istanton - anti-istanton separation increases). Ac- 
cording to the previous result on the essential spectrum, these (normalisable) modes, 
if there, would be embedded in the continuous spectrum. 

Non zero eigenvalues ±A for V correspond to a positive (doubly degenerate) eigen- 
value for the operator T)^. ^ The operator "D^ is substantially a Schrodinger oper- 
ator in R^, for which the problem of positive eigenvalues is well studied ([H2|). The 



following result (see [0 for the general case) is crucial for us: 

Theorem Let P = —D^D^ + c{x) (D^ = [d — igA)^) be a differential elliptic operator 
in euclidian R" with coefficents analytic almost everywhere, and such that the field 
strengths F^iy{x) = [D^,Di^] and c(x) go to zero at x 00 faster than j^, then P 



^It is a well known result (see [^) that for a self adjoint operator T, app{T") = io'pp(T)Y' 
{A"/A e cTppiT) (where app is the set of all eigenvalues of T). 
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cannot have positive eigenvalues. ^ 

Since the operator V'^ satisfies all the conditions required, the presence of positive 
eigenvalues can also be excluded. 



3. Large R/p 

In the absence of the fermion zero modes, the anomalous imaginary part in the 
forward elastic amplitude cannot be computed in a straightforward manner as in 
the instanton calculation of the anomalous production amplitude (|1.1|) , (|1.3|) . Still, 



at large R one expects physically that the generating functional should reduce to 
the product Z'^"'^ ■ Z^^\ where Z'^"'^ (Z'^^^ ) is the generating functional in the pure 
antiinstanton (instanton) background. 

Let us introduce complete sets of orthonormal modes {r]!;^^ and {C^*-*}, n = 
0, 1, 2, for the left-handed and right-handed fermions, respectively. They are eigen- 
states of and D^'W^'^ 0: 

D^''¥r:^ = fcmCi") (m = 0, 1, ...), = k^vl:^ (m = 1, 2, ...), 

D^Sl^^ = (m = 1, 2, ...), D«Cl:^ = Lr/» (m = 0, 1, ...), (3.1) 

where 

ko = lo = 0. (3.2) 
The functional integral can then be defined as: 

J V^Vij= n damdbn] (3.3) 

m,n=0 



^"The theorem of Ref.|47 refers to the case of Abehan backgrounds, but we beheve that the proof 
is general enough to be appHed to the non AbeHan case, too. 

The covariant derivatives D'^°'\ D^"^^ are defined with respect to the antiinstanton A^"-* and 
the instanton A^,^ of Eq.( |1.6| ). Accordingly the zero modes are those in the regular gauge (for the 
lefthanded mode) and in the singular gauge (for the righthanded one), respectively. See Appendix 
D. 

^^In Ref.[^ the valley approach has been applied also to fermions, to perform the functional 
integration. The authors also use a projector to select the anomalous intermediate states, which 
however is not explicity defined. It is not clear to us how such calculations can be done in practice, 
but we expect nontrivial corrections at finite R to the leading term ( 3.23| ). 



10 



X 



X) 



(3.4) 



m=0 



n=0 



As is clear from the way Eqs.( p.l| )-( p.4| ) are written, we first put the system in 
a large but finite box of linear size L ( such that L ^ R, p) so that all modes are 
discrete. After the derivation of Eq.( |3.9| ) below (i.e., after the sum over the complete 
sets is done), however, L can be sent to infinity without any difficulty. 

The two point function I{x, y) can be written as 

I{x,y) = det D {x\D-^\y) 

= detD{{x\a,0){a,0\D-^i,0){i,0\y) + ^ {x\a,m){a,m\D~^\i,0){i,0\y) 

+ '^{x\a,0){a,0\D~'^\i,n){i,n\y) + ^ {x\a,m){a,m\D~^\i,n){i,n\y)} : 

(3.5) 

the term proportional to the product of the zero modes has been singled out. We 
wish to compute Eg. (|3.5| ) at small p/R. To do this, first let us write the operator D 
in the above basis: 

f d Vi ... Vn . . . \ 



D 



wi Xn . 



■ J 



where we defined 



(3.6) 



d = (L))oo = (2,0|L'|a,0) = (i,0|C|a,0), 
Vn = {D)on = {'i,0\D\a,n) = {i,0\B\a,n), 
Wm = {i,m\D\a,0) = {i,m\C\a,0), 
X.mn = {i,m\D\a,n). 
Operators C = C^a^; B = B^a^ are defined by: 

Jjivalley) ^ jj{a) ^^^^ ^ 

The explicit expression of these quantities are found in Appendix D. 



(3.7) 
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The idea is that the matrix elements involving either of the zero modes, d, f„, 
are all small by some overlap integrals (see Footnote 13) while the matrix elements 
X are large because the wave functions of non zero modes are extended to all over 
the spacetime. 

The inverse matrix is given by: 

iD-\o = l/id-vX~^w) 

= + d^'^Vm{X'^)mnWn^ ; 

iD-%n = {X --w(E)v)-' = X-\l--w®vX-Y' 
d d 

= {X'^)mn + d~^{X'^)mlWlVk{X'^)kn H , 

(D-l)o„ = -d-\{D-')ln, 

{D-')mo = -{D-')ooXj,Wk, (3.8) 

where X~^ is the inverse of the submatrix X in the space orthogonal to the zero 
modes. 

Inserting Eq. (|3.8| ) into Eq. (|3.5| ) and after some algebra (see Appendix B for deriva- 
tion) one finds a remarkably simple (and still exact) expression for I{x,y): 

I{x,y) = detX{(x|a,0) - {x\X~^C\a,0)} {{i,0\y) - {i,0\BX-^\y)} 

+ detD {x\X-^\y). (3.9) 

Unfortunately, the inverse X~^ is not known, and our present knowledge does not 
allow us to compute it perturbatively in any small parameter. 

In Ref . , we used a perturbative formula for X^^ in terms of the propagator 
in the anti-instanton background S, to get an expansion for I{x,y). A closer look 
at various terms arising from such an expansion, however, has revealed several diffi- 
culties with this procedure. First, there are difficulties in the application of the LSZ 
procedure on some terms, caused by the infrared behaviour of S. Also, the expansion 
resulting from Eq.(12) of ||31[] (obtained from Eq. ( p.9| ) by use of the above mentioned 
formula for X~^ ) turns out not to be an expansion in p/R contrary to the incorrect 
claim made there. 



Nonetheless, we believe that Eq. (|3.9|) displays the main features of the two point 
function in the valley background correctly. The effect due to the zero modes is sepa- 
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rated explicitly and everything else is expressed by the smoother two point function, 

S',^^={x\X-'\y). (3.10) 

In order to find the amplitude, one must apply the LSZ reduction on the four 
point function I{x,y)I{u,v). This requires only the knowledge of S'^ y at large x and 

y- 

We assume S'^ y to behave at large x and y (with Xj and Xa fixed) as 

S'^^y ~ U\x)SAx,y)U{y), (3.11) 



U{x) 



c^yx Xa)^ 

\J{X- XaY 



where Sp is the free Feynman propagator. This behaviour is suggested by the fact 
that the valley field has a pure gauge form at large x, 

^(valley) _ 1 ^ _ Q^}^y (3_;L2) 

9 X 

To proceed further we assume that 

detX/det(9 = const. (3.13) 



Next the ratio detD/ detX can be estimated as follows (see (|3.8|) ): 

detZ) 



detX 

where use was made of 



{{D'^)oo)-^ ~ c/ ~ const. p'/R^ (3.14) 



d = Doo = {t, 0\C\a, 0) = / Cl;\zrC{z) r,t\z) ~ / R\ 



Combining Eq. (|3.13|) and Eq. ( |3.14|) gives 

detD 
det(9 



p^/R^. (3.15) 



With Eq. (|3.15|) and Eq. (|3.11|) in Eq. ( p.9|) one can extimate the amplitude and the 



leading contribution to its anomalous part. A naive application of the standard LSZ 
procedure, i.e., 

lim / ci1re^9-^'M(g) a^d^ {T^{x) ■ ■ ■) • (3.16) 

o^— >0 J 
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however leads to a number of difficulties. On the one hand, a large correction is found 
from the second term of (i,0|(l — BX~^)\y). The anti-instanton zero mode (in the 
regular gauge) {x\a,0), on the other hand, has not the right asymptotic form of the 
free propagator. 

All the problems however can be attributed to the inappropriate choice of the 
gauge, Eq. (|1.6| ). Indeed, we have already noticed that behaves asymptotically 

as a pure gauge field, Eq. (|3.12|) . It is then not surprising that the naive LSZ fails: 
the fermions never become free, whatever distance they travel, a situation somewhat 
reminiscent of the Coulomb scattering. 

Note that if the fermions were moving in a single instanton field, it would have 
been sufficient to work in the singular gauge, which also satisfy the Lorentz gauge 
condition. For the valley field instead, a gauge transformation to a sort of double 
singular gauge 



in which both instanton and antiinstanton do drop off faster at infinity, takes us out 
of the Lorentz gauge. 

We shall by-pass this problem by computing the Green function in the original 
gauge Eg. (|1.6|) and by transforming it to the more physical gauge Eq.( |3.17D only at 
the end of calculation, before applying the LSZ reduction. Note that the two point 
function I{x,y) = {Tip{x)il!{y)) transforms covariantly. 



valley) 



jj^ivalley) ^ Iq^^^U\ 



u 




(3.17) 



I{x,y) 



U{x)I{x,y)U\y). 



The correct LSZ procedure is thus: 





(3.18) 



for the initial fermions; 





(3.19) 
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for the final fermions. 

We analyse first the contribution 

{{x\a, 0) - {x\X-^C\a, 0)} {{t, 0\y) - {i, 0\BX-^\y)} (3.20) 

to the two point function (|3.9|). 

Considering the term {{x\a,0) — {x\X~^C\a,0)} and applying the LSZ amputation, 
we obtain: 

dU{x) 4%) - U{x)C{x)r]i%), (3.21) 

where we used ( p.ll| ) and d{SF)x,y = S{x — y). Note that the righthanded zero mode 
Tjo^Xx) is automatically transformed into the singular gauge form U{x) riQ"'\x) which 
has the correct asymptotic behaviour. 

Also, the large correction arising from the second term of {i, 0|(1 — BX^^)\y), due 
to the bad asymptotic behaviour of -B^ = _d(^^'^''^?') — Z) W gets automatically cancelled. 
To see how this occurs, first write 



Then 



b: = o{i/x') 



asymptotically. Then 

[{z,0\y)-{i,0\BX-'\y)]U^d 
~ {i,0\y)dU^ + {i,0\y){U^dU)U^ 
- {i,0\y)[{U^dU) + B']U^SFd 

= C%ydU\y) + C%yB'{y)U\y), (3.22) 

where use was made of the identity, {Sp)x,ydy = 6{x — y). 

Multiplying (|3.21| ) and ( p. 221) and doing the Fourier transform prescribed by LSZ 



we obtain the contribution of the first term of ( |3.9| ) to the amplitude. The leading 
term in this amplitude, containing exp(ip-a;a) exp(— 'ig-a;^), arises from the product 
of the zero modes in Eg. ( p.20D , as seen from the asymptotic behaviour of the latter. 
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This contribution to the amphtude must obviously be considered as anomalous. All 
remaining terms in ( p.20|) give a contribution of order 0{p/R) or less to the total 
amplitude and so are negligible for large R. In particular, second term of Eq. (|3.22|) 
gives a 0{p/R) contribution which is still proportional to exp{ip ■ Xa) exp{—iq ■ Xi): 
such a correction to the anomalous part of the elastic amplitude can become important 
when R p. 

The second term of (|3.9| ), with the assumption ( p.llD for S'^y, on the other hand, 
gives clearly a non anomalous contribution to the amplitude. Also eventual subleading 
asymptotic terms of S'^ y contributing to anomalous part would be suppressed by the 
small factor ( p.l4| ) compared to the leading anomalous term previously found. 

Recapitulating, the leading anomalous contribution ( for R/ p^ 1) is essentially 
given by 

J(™)(x,i/) ^ const. r/S"\x)cT(l/)*. (3.23) 



Inserted in the four point function Eq. (|1.9| ), it leads to the "anomalous part of the 
forward elastic amplitude" , required by unitarity and chiral anomaly. 

To reach the above conclusion really, however, we must make one further check. 



For the first term of Eq. (|3.9|) to represent the anomalous process, the intermediate 
state must contain fermions satisfying the instanton selection rule, Eq. ( |1.3| ). 

That this is indeed so can be seen from Eq. ( p.lSl) . The functional integration 



yields, for each flavour = 3, 4, ■ ■ ■ Np) one factor of 

det det 9 ~ pV^^ 

but this is precisely the factor expected for a left-handed fermion, produced at the 
instanton center (with amplitude p ), propagating backward to the antiinstanton (with 
amplitude, ~ ^-/{xa — XiY = 1/R^ ) and absorbed by the latter (with amplitude, p 
), see Fig. 1. (This provides an a posteriori justification of the assumption Eq.( |3.13|) 
for detX.) 

To be even more explicit, suppose that we are going to observe the final state by 
setting up an appropriate detector. This would correspond to introducing a source 
(or sink) term for each flavour, / dx J'^^{x)il)i{x) + /i.e., and taking the first derivative 
with respect to the sources. This would produce pairs of zero modes (as in Eq.( p.23| )): 



Clearly suppression is due to the appearence of product of functions with small overlapping 
support or to the explicit factor in (see appendix D for its explicit form). 
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the fermions required by the topological selection rule are indeed there, as long as 
R/ p is large (see Fig. 2). 

In this connection, it is amusing to note how a sort of dilemma between the 
factorisation in the flavour of the functional integration Eg .( 11.91) and the topological 
selection rule, is solved in the semi-classical approximation we are working with. 
Consider as an illustration the forward elastic amplitude with two initial particles 
of the same flavour, 1 + 1 — > 1 + 1. In this case no anomalous piece is expected 
to be present. To see this, let us note first that due to the Grassmanian nature of 
the fermion fields and sources one gets the amplitude with the Fermi-Dirac statistics 
automatically built in. Due to its factorized form, contribution of Eq. (|3.23| ) to the 
four point function cancels out by antisymmetrisation m. x ^ u\ y ^ v. 

As a further check, consider the case of massive fermions (we have so far re- 
garded them massless). For definiteness, consider Dirac type masses: we write 
ipi = ip^i^^ ; ip2 = = '^l^]4'a = etc. The contribution to Z of a sin- 

gle Dirac flavour is given by: 

Z = 

where the matrix acts between = [1^1,1^2) and ip = {ip2,'ipi)- 

At large R/pwe know that the first column and the first row of D and D, regarded 
as matrices in our basis (Eq. ( |3.1| )), become small due to the tiny overlap integrals, 
while the rest of the matrices remains large. It then follows that 

lim ^ = const, mm*, (3.24) 

R/p~*oo 

which is consistent with the clustering for massive but sourceless fermions. This 
argument further implies that Eq. (|3.24| ) is a good approximation when 




i?3 

4. Small R/p. 

4.1. Instanton -Anti-instanton Melting. 



< |m|. (3.25) 



The approximation p.23| ) fails at small R/ p < 1 for obvious reasons. 



17 



In particular, in the limit R/ p — > 0, the classical field A^^"'''''^^^ reduces to the 
trivial, perturbative vacuum. In the free theory, the procedure adopted above is still 
formally valid, but the appearance of the product of the zero modes (as in Eq.( p.5|) ) 
is of course a fake, the total two-point function being simply a^{x — y)^/{x — y)^. 
This means that at certain R/p ~ 0(1) the leading anomalous term Eq.( |3.23| ), must 
be effectively cancelled by contributions from other terms. Altough we do not know 
exactly at which value of R/p this occurs, there is a fairly good indication that such 
a transition takes place around R/p ~ 1. 

We have in fact computed numerically the integral up to the Euclidean time X4 
of the topological density, 

C(X4) = -£ll ^'^Ye^Tr F^^F^, = ^^csix,) ~^^csi-oo), (4.1) 

as a function of X4 for several values of R/p, for the valley background of Eq. (|1.6|) . 
(See Fig. 3.) J\fcs{x4,) is the Chern Simons number 

Afcsix,) ^-J d'x^e^'^'Ti {F.^Ak - ^AA^Ak). (4.2) 

The instanton and antiinstanton are situated at (0,i?/2) and at (0, — i?/2), respec- 
tively. 

In Fig. 4 we plot also the behaviour of the maximum of each curve, corresponding 
to C(0), as a function of R/p. The behaviour of C(0) is powerlike both at large and 
small R: 

C(0) ~ l-12{p/R)\ i?/p>l, 
C(0) ~ liR/pf, i?/p«l. 

Actually, the exact expression for C{x4) in terms of X4 can be found by using 



conformal transformations, |39|. In particular it can be proved that 



C(0) = 3(^)^-2(^)3 (4.3) 
z + 1 z + 1 

where z is defined in (11.7]). 

It can be seen from Fig. 3 that the topological structure is well separated and 
localised at the two instanton centers only at relatively large values of R/ p, R/ p > 10. 
At very large R, C{x4) approaches the product of two theta functions. Vice versa. 
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for small R/ p < 1 the gauge field is seen to collapse to some insignificant fluctuation 
around zero, not clearly distinguishable from ordinary perturbative ones. In this case 



there will be no level-crossing ||33|, ^ hence no chiral anomaly. 

Fig. 3 and Fig.4 indicate that the instanton and antiinstanton start to melt at 
around R/ p ^ 5 and go through the transition quickly, the center of the transition 
being at around R/ p = 1 . 

The behaviour of C{x4) at small R suggests that at R/p < 1 one can compute 
I{x,y) by simply using the standard perturbation theory, with ^(^^"'''^y) as perturba- 
tion. Our results have several interesting implications. 



4.2. Fermion Number Violation in the TeV Region Scatter- 
ing Processes. 



According to the method developed in Refs. p, ||, |^, [TJ, [T^, taking into account 



the contribution of the Higgs field in the action and using the saddle point approxi- 
mation in the integrations over the collective coordinates, one finds a relation among 
the cm. energy ^/s and the saddle point values of the parameters R and p. It was 



found in Ref.[|l^], in a simple toy- model calculation which uses the valley trajectory 
of Eq.( |1.6| ), that at an energy of the order of the sphaleron mass the total action (at 
the saddle point) vanishes and that the 't Hooft suppression factor disappears. 

But at that point, the corresponding saddle point values of the instanton parame- 
ters are found to he p = 0; R/ p = (hence = 0) : the unsuppressed cross section 
should simply correspond to a non- anomalous, perturbative cross section. This was 



pointed out also by some authors [17 



What happens is that the anomalous term of Eq.( |3.23| ) disappears at some R/p. 
Our numerical analysis of the Chern Simons number shown above suggests that this 
occurs quickly at around R/p = 1. And this turns out to be precisely where the 
valley action sharply drops to zero: see Fig.5 taken from Ref. |15[|. We must conclude 



that the toy-model considered there shows no reliable sign of anomalous cross section 
(associated with the production of large number (~ 1/a) of gauge and Higgs bosons) 
becoming large. 

Of course, there is no proof that all sorts of quantum corrections to the inst- 
anton-induced process (|1.1| ) are effectively described by a classical background such 
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as Eq. (|1.6|) . But this is another issue. The point here is that there are as yet no 
calculations anyway which show that the anomalous process becomes observable in 
high energy scatterings. 

More generally, the consideration of this paper suggests that, in order for fermion 
number violating cross sections to become observable at high energies, a new mech- 
anism must be found in which the background field governing the elastic amplitude 
Eg .( [1.21) does not effectively reduce to A'^J-I' = 0. Note that this is necessary, whether 
or not multiple- instanton type configurations become important. It is difficult to 
envisage such a novel mechanism, not accompanied by some finite fraction of the 't 
Hooft factor. 



4.3. Theories without Fermions. 

It is often stated that fermions are not essential, as dynamical effects of instanton- 
induced cross sections can well be studied in a theory without fermions. This is 
certainly true, but this does not mean that the consideration of this paper is irrelevant 
in such a case. 

Quite the contrary. The crucial factor, 

exp{ipj ■ Xa) exp{-iqj ■ Xi) {j = 1, 2) 

associated with the external particles, in the case of process with fermions, appears 
upon LSZ amputation applied to the product of the zero modes of Eg. ( p. 231) . In the 



case of external gauge bosons, the same factor emerges as a result of the " semiclassi- 
cal" approximation, 

< TA^{x)A^{u)A^{y)A^{v) > 



A\^^\x)Ai^\u)A^;\y)A^\v) + .... (4.4) 



It is however clear that this approximation, which is probably good at large R/p, 
should fail at R/p < 1 just as Eg.(|3.23|) does. The first term of Eg.(|4.4|), the anoma- 



lous term, disappears precisely when the cross section is claimed to become observable. 

Thus the use of the theories without fermions does not change our conclusion: 
there is no evidence for the anomalous process to become observable in the TeV 
region scattering. The scarcity of our knowledge is actually even worse. In fact. 
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in the existing literature the substitution Eq. (|4.4|) is just made by hand : no proof 
of Eq. ([4.4| ) seems to be known. It is an interesting question how our approach (of 
Section 3) can be generahsed to the case of external bosons. 



5. An Apparent Paradox and its Resolution. 

It is quite remarkable that the fermion number violating term (Eq.( |3.23D ) at large 
instanton-antiinstanton separation, emerges without there being a single, dominant 
mode of the Dirac operators, D, or D. The result Eq.f ^.lSD for det D neither implies 



the existence of a particular eigenmode with eigenvalue, ~ p^/B?, nor requires that 
such a mode dominate over others. 

Indeed, for any finite R we have established in Section 2 the following. First of 
all, neither exact zero mode nor quasi-zero modes exist in the valley background. 
Secondly, there are many non-zero modes, definitely lying below p^/i?^, and forming 
a continuous spectrum, reaching down to 0. In particular, putting the system in a 
large box of fixed size L such that 

L:^R,p, (5.1) 

for the lowest lying modes with kR <^ 1 both the eigenvalue and the wave function 
differ little from the free spectrum, as shown perturbatively in Appendix C. 

(In passing, this shows that a single i — a pair in itself cannot lead to chiral sym- 
metry breaking, < i/jip >^ 0, in QCD: the latter necessarily requires an accumulation 
of eigenvalues towards (see Ref. [0). In the context of instanton physics, that 
would require something like the "instanton liquid" (Ref. [0).) 

The results such as Eq. (|3.23| ) and Eq. (|3.24| ) are thus collective effects in which 



many modes contribute together; no single mode of D or D plays any particular 
role. We then seem to face a somewhat paradoxical situation. At large instanton 
- antiinstanton separation, physics must factorise and we do find results consistent 
with such intuition. This is fine. The problem is that mathematics to achieve this 
looks very different from that of the usual instanton physics where a single fermion 
zero mode plays a special role. In our case, there is no hint even of the presence of a 
quasi zero mode. What is going on? 
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The key to the resolution of this apparent paradox is the inequahty, (|5.1| ). In 
order to be able to compute the S-matrix elements from the four point function (LSZ 
procedure), we are in fact forced to work in a spacetime region whose linear size is 
much larger than the physical parameters R and p, independently of the ratio R/ p. 
All fields must be normalised in such a box and functional integration defined there. 
This makes our system with instanton-antiinstanton background always distinct and 
not continuously connected to the system with a single (say) antiinstanton. The 
gauge field topology remains firmly in the trivial sector. 

In spite of our use of the "zero modes" r]^\ (q^ as a convenient device for calcu- 
lation, they are not a good approximate wave function for the lowest modes in the 
valley background, however large R/p may be. In fact, if one insisted upon using 
tiq""^ as the "unperturbed" state, one would discover that the effect of C is always 
nonperturbative and large (near Xi). ^ 

All this is clearly to be distinguished from the situation where we ask e.g., what 
the effects of an instanton on a distant planet are. ^ Such a case would correspond 
to the inequality opposite to (|5.1|) , 

if we restrict ourselves to the instanton-antiinstanton case. As long as we are inter- 
ested in physics inside our laboratory (or on the earth, anyway) both physics and 
mathematics are described by just ignoring the distant instanton, to a good approxi- 
mation. (Particles are produced and detected in the region of volume, -C -R^, the 
gauge field having the winding number —1, fields and functional integrations defined 
in the same "box", etc.) The effect of the distant instanton is a true, and negligibly 
small, perturbation in this case. 



We did verify that the standard perturbation theory apphed to H = DD = + H' 

with H' = D'-'^W + CD^"''^ + CC, yields Ai; = to aU orders, reflecting the topological stability of 
the fermion zero mode. (Another way to see this is to notice a supersymmetric structure underlying 
the system. We thank C.Imbimbo for pointing this out to us.) Such a result, however, is false in 
the case of topology- changing modification of the gauge field, in view of the result of Section 2. 
^^A somewhat related consideration is found in Ref. [P5|. 
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6. Conclusion. 



In this paper we study the fermions propagating in the background of instanton- 
anti-instanton valley. The spectrum of the Dirac operator is studied first, proving 
the absence of the zero mode for all values of the collective coordinates. We then 
study the fermion Green function at large R/ p (the instanton and anti-instanton well 
separated). This solves the problem raised in the introduction of unitarity versus 
chiral anomaly, yielding the leading anomalous part of the forward elastic amplitude. 
For small R our study remains somewhat indirect, being limited for the moment to 
the numerical study of Chern Simons number as a function of the time and of the 
ratio R/p. Nonetheless the latter has several interesting implications as was discussed 
in Section 3. 

Applied to the question of fermion number violation in the TeV region scattering in 
the standard electroweak theory, our results imply that there is as yet no theoretical 
evidence that such a process becomes observable, contrary to the claim of several 
papers. 

In the application to Quantum Chromodymamics, we hope that our study serves 
as the starting point of an improved study of the instanton liquid model for the QCD 
vacuum. 

Acknowledgment. The authors are grateful to D. Amati, M. Bertero, C. Imbimbo, 
H. Leutwyler, G. Morchio, M. Maggiore and V. Zakharov for interesting discussions. 



Appendix A: Essential Spectrum of the Dirac Operator in 
the Valley 

In this Appendix we choose for the valley field the 'clever gauge' introduced in 
Section 2: 

Aa = (Tp{af,au - O^ujC^T — 5- 

g 

^•^ = 7 \2 I 2 ~ 7 \r~r^ {v = x-Xi + y) 

[X — Xa) + P \^ ~ ^i) + P 

which is obtained from Eq. (|1.6|) by the gauge transformation, U = ^4^. is a 



bounded function and A^{x) ~ as x ^ 00, (which implies the uniformity in 
direction of x of convergency at infinity). The following theorem can then be used: 
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Theorem If A^{x) is a hounded function going to at oo, uniformly in x, then the 
essential spectrum of V = i'^^{d^ — igA^) in euclidian space is aess = (~oo, oo). 

Proof: We treat the vector potential term as a perturbation of the free Dirac opera- 
tor, and we prove that complete and the free Dirac operator have the same essential 



spectrum. It is a well known result |44, cap. IX] that the free Dirac operator 

(that is acting in the space of square- integrable four-spinors) , with domain -D(Po) = 
if^(R^) 13 is self-adjoint (and so closed) with spectrum 

a(I?o) = cress(^^o) = ("OO, +CX)) 

Thanks to the boundedness of vector potential the operator A : L^(R^)^ — > L^(R^)^ 
which acts multiplying spinors by g'y^Afj_{x), is bounded and self-adjoint, so that 
the operator V = Vq + A is self-adjoint (with D(T>) = D{T>q)). 

To prove the equality of the essential spectra we use the following theorem (see 



H, p.ll4 and p.ll6]) with T = and 1/ = A: 



Theorem Let T be a self-adjoint operator in a Hilbert space, and V a bounded self- 
adjoint operator. Also let V T'^-compact for some integer m > 0. Then 

aess{T + V) = (TessiT) 

An operator V is called compact relatively to a self-adjont operator T (shortly 
T-compact) if D{T) C D{V) and the operator V{T — fi)~^ is compact, for some 
fi ^ <y{T). We shall not give the definition of compact operator here (see |^), but 
just use the result that in R" an integral operator with kernel K{x, y) such that 

d'^x d''y\K{x,y)\^ < oo (A.l) 

is compact. 

In order to use this theorem with T = Vq, V = A, m = 4 and fi = —1 we must 
prove that the operator A{d^ + 1)~^ is compact. To do this consider the sequence 

^^We indicate with i?™(R") the space (called Sobolev space, see Q) of all functions which are 
in L^(R") together with their first m (distributional) derivatives. These spaces, if endowed with an 
adequate topology are Hilbert spaces. 

"'^^That is II Alp \\< B \\ \\ Vip, property which is equivalent to continuity. 

i^Note that, we use T = Va, therefore = 0". The following proof is inspired to |4|, p. 117]. 
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of truncated vectors potentials 

Anf,{x) = Af,{x) 9{n-\x\), n = l,2, ... 

which, due to the property of (uniform) convergence to zero at infinity of the vector 
potential, converges to uniformly in all x. Q This guarantees that also the 

sequence of associated multiplicative operators An = (77^A„^(a;) converges to A in 
the (bounded) operators space, and, by continuity of multiplication, the sequence 
Anid'^ + l)""*^ converges to A{d'^ + 1)~^ also. 

Because the space of compact operators is a closed space, it then suffices to prove 
that for all n the operators An{d^ + 1)~^ are compact. But these operators are integral 
operators with kernel: 

Kn{x,y) = A^{x) 9{n - \x\) G{x,y) 

where G{x, y) = S p^qi^y is the euclidian 4-dimensional Green function of 9^ + 1. 
It easy to check that the kernels Kn all satisfy condition ( |A.1| ), so that An{d^ + 1)~^ 
and A{d'^ + 1)~^ are compact. It means that A is Pg-compact and therefore we can 
apply the previous theorem to T) = Vq + A, obtaining the result: 

a{V) = aessiP) = (-00, +oo). 

We conclude that if T> has eigenvalues (with normalisable wave functions) they are 
embedded into continuous spectrum and they are not isolated. Q.E.D. 



Appendix B: Derivation of Eq. (|3.b|) . 



In the basis ( |3.1| ) write: 



D = 

Imposing DD^^ = 1 one has (summed indices are left implicit): 

dd' + vw' = 1 



d 


Vn 


Wm 


-^mn _ 


' d' 






VI 

nm 



(B.l) 



"'^^I.e. in the || ||oo norm. 

^"For operators that acts multiphng by a function, the norm in the operator space coincides with 
the II lino norm of that function. 
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wd' + Xw' = (B.2) 

dv' + vX' = (B.3) 

w0v' + XX' = 1. (B.4) 
Solving these equations for the primed quantities one derives Eq.( p.8| ). 
Also DD^^ = 1 gives: 

d'd + v'w = l (B.5) 

w'd + X'w = Q (B.6) 

d'v + v'X = (B.7) 

w' ®v + X'X = 1. (B.8) 
Using Eq. ( |B.8| ) and then Eq. ( [B.2| ) one finds: 

X'-i - = -w' ® {vX~^) = d'iX'^w) ® {vX"^). (B.9) 

Note also the relations 

w' = -d'X-^w, (B.IO) 

v' = -d'vX-\ (B.ll) 

following respectively from Eq.( |B.2|) and from Eq. (|B.7| ). Using Eq.( |B.9| - |B.ll| ) one 
obtains: 

{x\D'^\y) = {x\a,0)d'{i,0\y) + J2{^\0'^^Wm{hO\y) 

+ ^(a;|a,0X(i,r2|y) + ^ {x\a,m)X'^^{i,n\y) 

nj^O m,n^O 

= {x\a,0)d'{i,0\y) - d' {x\a,m){X'^w)m{i,0\y) 

m/O 



- d'J2{x\a,0){vX ^)n{i,n\y)+ ^ {x\a,m){X' - XjJ{i,n\y) 
+ {^\(^^^)^mn{hn\y)} 

= d'{{x\a, 0) - {x\X-'C\a, 0))((z, 0\y) - {z, 0\BX-'\y)) + {x\X-'\y) 

(B.12) 
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The inverse of a matrix M is M ^ = det M ^Cof M*, so that 

detDrf' = detX. (B.13) 
Substituting (|B.13|) in (|B.12D , one derives the expression Eq. (|3^ ) for 

I{x,y) = det D{x\D-^\y). 

Appendix C: Perturbative Approach to the Spectrum of 
DD. 

We consider the hermitean 'Hamiltonian' H = —DD acting in the space of func- 
tions normahzable in a box of size L ^ p,R (with periodic boundary conditions). We 
wish to study the spectrum of this operator with the ordinary perturbation theory 
of discrete spectrum (in the hmit of L — oo the spectrum of H is continuous, see 
Appendix C). 

For unperturbed hamiltonian we take Hq = —d"^, with eigenfunctions 

{x\k) = ^uik) e^'^^ A;=™ 

associated with eigenvalues E^'^ = k"^. With this choice, the perturbation is 

V = Ad + dA~A^ 

where A = —iga^A^^"'^^'^y\ We take the valley background in the 'clever' gauge pre- 
sented in Section 2, in which A is everywhere regular and behaves at infinity as 

A{x) ~ 4 
x'' 

(we neglect all spin and isospin factors that are not relevant to our approximate 
analysis) . 

Our aim is to show that the spectrum of H, for kR <^ 1 is very close to that of 
Hq, i.e. that perturbative corrections are small compared to E^^ and perturbation 
theory is well posed. 

Consider a generic unperturbed level E''^\ which is degenerate (apart the case 
/c = 0), thus first correction Ej^^ is one of eigenvalues of the matrix 

Vkk' = {k\V\k') (C.2) 

~ t{k\q^A^ + 2k^A^\k'), (C.3) 
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where we have neglected the second order term and we have set q = k — k'. 
For our purposes it is sufficient to consider 

m&x{k,q) 

Vkk' = V{q) ^ /(g) 

where we defined I{q) = dhce^'^^ A[x) (from here on we consider only moduli of 
vectors). Note that for degenerate levels fc^ = fc'^ it follows that < |g| < 2\k\. 

For values of q such that gL ~ 1, the main contribution in the integration comes 
from the region at infinity, where we can use Eq. (|C.l|) for A, arriving at the estimate 

7(0) = " ~ 0^' " = - 

(we used here a spherical box). 

The Riemann-Lebesgue theorem guarantees also that, for a fixed L, 

im / 



lim / rfW^A(x) = 0, 

q^oo J 

which implies that 



lim ^ = 0. (C.5) 



Equations (|C.4|) and (|C.5|) lead to the relation 



I ^(g) I ^ max(A;, q) , . _ , n ^ 

which we assume to be valid for all q. 

The degeneracy of each level /c = ^ is of order M ~ (ra = |n|), thus, using 
Eq.( P.6|) , we can roughly bound the eigenvalues of perturbation matrix Vkk' by: 

l<l<M2|F(0)|~^, j = l,---M. (C.7) 

For the perturbation theory to be good, the first order correction must be small 
compared both to the zeroth order value and to the gap with the nearest nondegen- 
erate energy level: 



Using the bound Eq.( |C.7| ) as an extimate for l-E'^J'*!, we see that these conditions are 
indeed satisfied for <^ ^, i.e. substantially for kR ^ 1, as we claimed. 
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Also, from Eq.(^) it follows that \Vk'k\ < \k'^ - k"^\ for all k', 

are eigenfunctions of matrix Vk'k and of course differ from unperturbed eigen- 
functions, due to degeneracy), substantiating our claim about smallness of corrections 
to the wave function. 



Appendix D: Conventions and Formulae 

In this Appendix are collected explicit expressions of many quantities frequently 
used in the text. 

We use the definition of covariant derivative: 

= (9 — igA'^)^; c = valley, i, a 
where external fields are given by: 

A (valley) _ /l(a) i A{i) i Aiint.) 



g [X - Xa) + p 

if- r N {X - Xi)yP 



g [x: ~ XiYiXx - XiY ^ p^) 

Aiint.) ( ^ ^ X \\{^~^i~^y)f {x — Xi)y^ (T\ ^\ 



and 



y = -Rl{z-l)- z = {R' + 2p' + jR^ + 4p^R^)/2p'; R^' = {xi - Xa) 



Note that, for > 1 \y\ ~ ^. 

In the text we use the definitions: 



D = a^Dl^^^^y-, D = a^Dl^^^^y 
D' = a,Dl- D- = a^Dl- c = i,a (D.2) 
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where we adopted the convention 



In spite of the indistinguished use of the matrices o"^, cr^, both for the colour (as in 
(|D . 1| ) ) and spin (as in (p.2| )), the reader will easily distinguish them from the context. 



The fields C^, are defined by 

^valley ^ jja ^ ^ ^ ^ JJ^^ ^ ^ ^ 

SO that 

Expressions of zero modes of operators D^"-* and D*^*) are: 

IT 

^ • I {x — XiY 



where we defined 



(x — Xc)^ + 

for c = i,a. 

Asymptotic behaviour of various quantities is: 



SUV / / 



[x — XiY R{x — Xi)^ 



However 



with f/a = ""'^(^ ^"^^ which shows that the dominant term of B^ at infinity has a pure 
gauge form. 
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Figure Captions 

Fig. 1. Each flavour j = 3, 4, • • • Np contributes a factor p^/R^. 
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Fig. 2. The propagators for j = 3, 4, ■ ■ ■ Np, cut at the intermediate state. 

Fig. 3. C{x4) versus x^/ p for = 10 (outmost curve), ^ = 5, = 2 (middle), 
= 1 and = 0.5 (innermost curve). 

Fig. 4. C(0) as a function of ^. 

Fig. 5. Valley action -^S'""-'-'-'^^ as a function of - (taken from |T5| ). 
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